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Computational fluid dynamics (CFD) analysis has been performed to study the plume 
effects on sonic boom signature for isolated nozzle configurations. The objectives of these 
analyses were to provide comparison to past work using modern CFD analysis tools, to 
investigate the differences of high aspect ratio nozzles to circular (axisymmetric) nozzles, 
and to report the effects of underexpanded nozzle operation on boom signature. CFD 
analysis was used to address the plume effects on sonic boom signature from a baseline 
exhaust nozzle. Near-field pressure signatures were collected for nozzle pressure ratios 
(NPRs) between 6 and 10. A computer code was used to extrapolate these signatures to a 
ground-observed sonic boom N-wave. Trends show that there is a reduction in sonic boom 
N-wave signature as NPR is increased from 6 to 10. The performance curve for this 
supersonic nozzle is flat, so there is not a significant loss in thrust coefficient as the NPR is 
increased. As a result, this benefit could be realized without significant loss of performance. 
Analyses were also collected for a high aspect ratio nozzle based on the baseline design for 
comparison. Pressure signatures were collected for nozzle pressure ratios from 8 to 12. 
Signatures were nearly twice as strong for the two-dimensional case, and trends also show a 
reduction in sonic boom signature as NPR is increased from 8 to 12. As low boom designs 
are developed and improved, there will be a need for understanding the interaction between 
the aircraft boat tail shocks and the exhaust nozzle plume. These CFD analyses will provide 
a baseline study for future analysis efforts. 


Nomenclature 


Cf g 

computed thrust (from CFD)/ideal thrust 

NPR = 

nozzle pressure ratio = PJP 0 0 

P 

local static pressure, psia 

Pt 

total pressure in nozzle 

Poo 

free-stream static pressure 

AP/Poo = 

(P-PJ/Pao 

A P 

P-Poo 

x = 

axial distance, in. 

D 

test nozzle diameter, in. 

x/D 

nondimensional axial distance from jet simulator nose cone 

I. Introduction 


S upersonic aircraft have operational restrictions over populated areas because of the disturbance caused by the 
sonic boom. Research is needed to enable the reduction of sonic boom for supersonic cruise over land. The sonic 
boom is caused by the coalescence of shock waves formed off the aircraft. Recent work has been performed to 
demonstrate reduction of the sonic boom N-wave for supersonic aircraft. Examples are the Supersonic Boom 
Demonstrator 1 and the Quiet Spike 2 program. These studies reduce the magnitude of the N-wave at the nose of the 
airplane by shaping the aircraft. A sample N-wave sonic boom signature is shown in Fig. 1. Work is also needed to 
reduce the magnitude of the N-wave generated by the aft airplane components, including effects of the exhaust 
nozzle plume. 
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Fig. 1. Sample N-wave sonic boom signature. World’s first ground-level 
shaped sonic boom recording in blue, baseline F-5E N-wave in red. 
Shaped sonic boom demonstration (SSBD) flight 9, August 27, 2003 (from 
ref. 1). 


Past work has determined that the sonic boom N-wave signature is based on aircraft volume and lift. 2 3 Also, 
previous studies examine the effect of various components, such as engine nacelles. 4 More recent work has 
determined that the sonic boom signature can be extrapolated from the near-field pressures of an aircraft as 
measured near the vehicle. 5 An analysis tool called “PCBoom” was developed using these theories to extrapolate 
pressures near the aircraft to an estimated ground-observed sonic boom. 6 

Previous work by Barger and Melson 7 show that highly underexpanded nozzle plumes have a large effect on the 
sonic boom signature. Also in previous work by Putnam and Capone, 8 a variety of nozzles were tested from a fully 
expanded mach 1.7 nozzle to a fully expanded mach 2.9 nozzle. Their study was conducted in a wind tunnel, where 
near- field pressure measurements were made at 1 nozzle diameter away from the test nozzle. Wind tunnel conditions 
were mach 2.2 at a simulated altitude of 50,000 ft. Inspection of their data demonstrated two trends, (1) near- field 
pressures are reduced when the nozzle is operated near its design point and (2) near-field pressures are reduced 
when the nozzle design point was matched to the flight conditions. Based on the trends from Putnam and Capone, a 
baseline nozzle was selected for a two-dimensional computational fluid dynamics (CFD) study. The nozzle 
configuration was selected to be operated close to the design point. For this study, “Nozzle 6” was selected as the 
baseline mach 2.0 nozzle in a simulated mach 2.2 flow field. 

The objectives of this analysis were to provide comparison to past work using modem CFD analysis tools, to 
investigate the differences between high aspect ratio nozzles and circular (axisymmetric) nozzles, and to report the 
effects of underexpanded nozzle operation on boom signature. 

Modern CFD tools were applied to the same nozzle configurations as tested by Putnam and Capone to provide a 
baseline for future nozzle plume studies. The CFD solutions were extrapolated to the far field to provide a baseline 
sonic boom N-wave. Additional CFD solutions were performed on a sample high aspect ratio nozzle configuration, 
using the baseline for comparison. 

During the analysis, it was determined that overexpanded and underexpanded operation of the nozzle have an 
effect on the N-wave boom signature. This trend was consistent for both the axisymmetric baseline and the high 
aspect ratio nozzle. This report summarizes the effect of an underexpanded nozzle plume on the N-wave boom 
signature for an isolated nozzle configuration, and demonstrates the feasibility of reducing the magnitude of the 
sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock stmcture. 
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II. Analysis 


A. Computational Domain and Grid 

1. Baseline nozzle 

The baseline nozzle chosen from the experiment conducted by Putnam and Capone is shown in Fig. 2, and the 
wind tunnel model used for testing is in Fig. 3. The computational domain for the baseline nozzle, Nozzle 6, is 
shown in Figs. 4 and 5. The domain consists of the axisymmetric nozzle with geometry from NASA TN-D-5553, 
and had an external boat tail angle of 5 degrees; nozzle geometry is shown in Fig. 2. The computational domain 
extends 18 nozzle diameters downstream of the nozzle exit, and 1.48 nozzle diameters above the nozzle outer wall. 
The length of the computational domain was extended beyond the experimental measurements in an attempt to 
capture additional effects and plume pressure signatures downstream of the nozzle exit. This grid was created with 
“Gridgen” software, and the first gridpoint off the nozzle wall was placed at wall coordinate distance, y + , of one. 
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Fig. 2. Nozzle configuration from NASA TN-D-5553. 
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The grid for Nozzle 6 is detailed in Figs. 4 and 5, and in Table 1. The overall grid was designed to be dense, 
especially in the free-stream zones, to capture the propagation of shocks off the nozzle geometry and plume. A grid 
density study validated the grid requirements, and is described later in this report. Based on the grid study, the 
baseline nozzle analysis used a two-dimensional grid of 51 1,299 gridpoints. 


Table 1. Baseline Nozzle 6 grid summary. 


Zone 

Name 

i 

J 

Gridpoints 

1 

Nozzle Inlet 

121 

161 

19481 

2 

Nozzle Exit 

419 

161 

67459 

3 

Free stream 1 

207 

103 

21321 

4 

Free stream 2 

111 

293 

32523 

5 

Free stream 3 

201 

293 

58893 

6 

Free stream 4 

311 

333 

103563 

7 

Free stream 5 

311 

333 

103563 

8 

Free stream 6 

311 

336 

104496 


Total 

511299 


2. High aspect ratio nozzle 

A high aspect ratio slot nozzle configuration was modeled to investigate the utility of high aspect ratio nozzles 
for reduction of sonic boom signature. The grid was modified from the axisymmetric grid used for Nozzle 6 with the 
same grid density. The internal nozzle contours were modified using a method of characteristics code to design a 
mach 2.2 slot nozzle, with an aspect ratio of 8:1. To simplify the analysis, the CFD simulation was performed on a 
two-dimensional slice through this nozzle. Geometry for this nozzle was a simple expansion ramp design with a 
2-in. throat height, a 3.88-in. exit height, and an overall length of 18.02 in. 

B. Computational Algorithm 

7. Wind-US 

Wind-US is a computational platform that numerically solves various sets of equations governing physical 
phenomena. 9 The code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along 
with supporting equation sets governing turbulent and chemically reacting flows. The current version of Wind-US is 
2.0. Wind-US was used with the modified second-order Roe upwind scheme for stretched grids, implicit time 
stepping with a Courant-Friedrichs-Levy (CFL) number of 1.0, and the shear stress transport (SST) turbulence 
model. No grid sequencing was used. Boundary layer thickness was specified on the nozzle outer wall in the mach 
2.2 free stream flow to match the boundary layer as measured in NASA TN-D-5553. 

C. Results 

7. Baseline nozzle 

Results for near-field pressures, which will be used to compute sonic boom signature, are typically plotted as 
AP/Pao. This is the primary input for the PCBoom code. Two computations of AT^/TU, were performed for this study. 

Both the current CFD results and the NASA TN-D-5553 wind tunnel results were performed at mach 2.2 and a 
simulated altitude of 50,000 ft. In the wind tunnel, near- field pressure measurements were made at 1 nozzle diameter 
away from the test nozzle. Figure 3 shows the wind tunnel model, where an x/D of zero is the front of the nose cone. 
The first effort was to validate the CFD. Data was taken from NASA TN-D-5553 Nozzle 6, which was digitized 
from the report and displayed in Fig. 6. A CFD data set was obtained at a nozzle pressure ratio (NPR) of 8, which 
would provide data near the design point of this nozzle, and would match the experimental test conditions. Figure 7 
shows contour plots of AP/P ^ and mach number for the CFD simulation. For a comparison of CFD results to 
experimental results, a cut was made through the CFD grid, also at 1 nozzle diameter away from the test nozzle as 
was done in the experiment. Cut planes are also shown in Fig. 7. 

The AP/^ pressure signatures were compared to the experimental data. The full results from the CFD are shown 
in Fig. 8, and a subset for direct comparison is provided in Fig. 6. The CFD computed the first AP/FC underpressure 
of -0.21 at an x/D of 8.55, followed by an overpressure of 0.11 at x/D of 8.7. This can be compared to the 
experimental data where the first AP/P ^ underpressure was measured as -0.2 at x/D of 8.6 and an overpressure of 
0.1 at x/D of 8.7. The difference in the maximum and the minimum pressures are 10 and 6 percent. Differences are 
attributed to the difficulty in surveying pressures in a supersonic wind tunnel. Problems can include effects from the 
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wind tunnel model supports and the ambient pressure changes in the test section . 10 Expansion waves from the wind 
tunnel or strut support can affect the initial expansion wave at the nozzle boat tail. CFD predicts a stronger 
expansion, which results in a stronger shock. 



Fig. 6. Near-field pressure profile AP/P ^ at 1 nozzle diameter from 
baseline Nozzle 6 axial positions x/D from 6 to 10 comparisons to NASA 
TN-D-5553. 




Fig. 7. (a) Baseline Nozzle 6 A P/P^ for the design NPR 8 at mach 2.2. (b) Closeup of baseline Nozzle 6 AP/P ^ 
for the design NPR 8 at mach 2.2. (c) Closeup of baseline Nozzle 6 mach contours for NPR = 8 and mach 2.2. 
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Fig. 8. Near-field pressure profile AP/P^ at 1 nozzle diameter from baseline Nozzle 6. 
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Fig. 9. Far-field pressure profiles A P/P^ at j = 150 in. from nozzle centerline. 
Baseline Nozzle 6 NPR 6 to 10. 


Figure 8 shows that beyond an x/D of 10.2, there was a second overpressure predicted by CFD analysis, which 
was not captured in any of the testing performed for NASA TN-D-5553. During testing, pressure surveys stopped at 
an x/D of 10. 

Data was taken from the CFD solution to use in the PCBoom code to extrapolate to an estimated ground sonic 
boom pressure signature. PCBoom requires the user to provide a pressure signature far away from the aircraft or 
nozzle to consider the source as a point source. This need is balanced with the need to keep a CFD analysis grid to a 
reasonable size and to obtain adequate propagation of external shock waves. Some trial and error was required on 
grid size to obtain reasonable results; ground boom signatures extrapolated from a location that was too close to the 
nozzle produced erratic results. 

A cut was made through the data at 150 in. from the nozzle centerline, or at 9.8 nozzle diameters, and this data is 
in Fig. 9. A series of nozzle pressure ratios were run, from 6 to 10, to develop a performance curve for this nozzle, 
and the AP/P X for each of these conditions is plotted. From this plot, the trend shows peak AP/P underpressure is 
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reduced as the NPR increases from 6 to 10. Increasing NPR increases the underexpansion in nozzle flow and causes 
the plume to expand, reducing the boat tail expansion wave and recompression shock, illustrated in Fig. 10. 

The data at 150 in. was used as an input to a subroutine of PCBoom, called “foboomv.” Flight conditions were 
50,000 ft at mach 2.2 for straight and level flight through a uniform atmosphere. The output is the sonic boom 
signature for all angles under the aircraft. This study used the signature at zero degrees, directly underneath the 
aircraft. In this case, there is no actual aircraft signature, so the signature is just from the nozzle and nozzle plume in 
simulated flight. Figure 11 shows the estimated sonic boom N-wave as perceived at ground level. The trend is 
similar from the pressure signatures in Fig. 9. The peak-to-peak strength of the computed sonic boom N-wave is 
reduced with increasing NPR. 


NPR 6 



1.2 1.599 1.998 2.397 2.4 

Mach number 


Expansion 



NPR 10 



1.2 1.5 1.8 2.1 2.4 

Mach number 



Slightly overexpanded Slightly underexpanded 

Fig. 10. Effect of nozzle operating conditions on the local shock waves. Baseline Nozzle 6 at NPR 6 and 10. 



Fig. 11. Estimated sonic boom signature at ground through an ideal 
atmosphere. Baseline Nozzle 6 at NPR 6 to 10. 


8 

American Institute of Aeronautics and Astronautics 


The sonic boom signature for the design point at an NPR of 8 differs in the forward part of the N-wave boom 
signature. The data input and processing with PCBoom was the same for all cases, and the difference does not affect 
the overall trend. 

2. High aspect ratio nozzle 

Following the baseline analysis, a simplified two-dimensional CFD analysis was performed to demonstrate if 
high aspect ratio nozzle configurations would provide for a favorable distribution of near-field pressures in the 
nozzle and plume system, which would reduce sonic boom. If favorable, a high aspect ratio nozzle would enable a 
highly integrated configuration, as long slot nozzles could blend into a wing. Potential advantages include a low- 
profile integration of the exhaust nozzle into a low boom vehicle concept. As a result, comparisons were done of a 
high aspect ratio nozzle to the baseline axisymmetric nozzle configuration. The high aspect ratio nozzle was 
designed using a method of characteristics code to match the exit area of the baseline nozzle. A high aspect area 
ratio was selected at 8:1. More aspect ratios would be analyzed if a favorable near-field pressure distribution was 
computed. 

A two-dimensional nozzle was designed using a method of characteristics computational code. The nozzle was 
designed to have the same area as Nozzle 6 from the previous axisymmetric study, and an aspect ratio of 8:1 was 
chosen. The nozzle was designed for an NPR of 8 at mach 2.2, and the CFD simulation was simply a two- 
dimensional simulation to study the far-field pressure and shock patterns. No end effects or three-dimensional 
effects were modeled for this study. 

As before, a cut was made through the data at 150 in. from the nozzle centerline, and this data is in Fig. 12. A 
series of nozzle pressure ratios were run, from 8 to 12, and the dP/P for each of these conditions is plotted. The 
shock patterns from the two-dimensional nozzle were much stronger than the axisymmetric case. Figures 13 and 14 
show the difference in expansion waves which are caused by the nozzle boat tail and the plume. As NPR is 
increased from 8 to 1 1 , the plume expansion shock interacts with the boat tail expansion shock. This can be seen in 
the AP/Pao signatures in Fig. 12. As NPR is increased from 8 to 12, the underpressure is changed from a negative 
value of -0.085 to -0.03. 



x , in. 

Fig. 12. Far-field pressure profiles AP/P ^ at y = 150 in. from nozzle 
centerline. Nozzle 6 slot nozzle NPR 8 to 12. 
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Fig. 13. (a) AP/Poo for Nozzle 6 slot nozzle configuration, design NPR 8. (b) Closeup of AP/P ^ for Nozzle 6 slot 
nozzle configuration, NPR 8. (c) Closeup of mach number for Nozzle 6 slot nozzle configuration, NPR 8. 




Fig. 14. (a) AP/Poo for Nozzle 6 slot nozzle configuration, NPR 11. (b) Closeup of AP/P ^ for Nozzle 6 slot 
nozzle configuration, NPR 11. (c) Closeup of mach number for Nozzle 6 slot nozzle configuration, NPR 11. 
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These pressure signatures show very strong AP/P signatures created by the two-dimensional geometry, and are 
nearly double that of the axisymmetric case. As a result, the PCBoom ground pressure signatures are not reported 
here, as the ground pressures were unreasonably high. A full three-dimensional simulation of this nozzle may 
produce more realistic ground pressure signatures. However, it seems reasonable to assess that the local expansion 
shocks generated by two-dimensional nozzles are locally stronger than that of an axisymmetric case. This result is 
consistent with the differences between shock structures created by a wedge and a cone. For low boom design, 
stronger shocks are undesirable. However, these results may still be beneficial in the study of highly integrated low 
boom aircraft and propulsion systems. Also, the two-dimensional solution used for this study may not accurately 
capture the three-dimensional effects of a high aspect ratio nozzle. Future small-scale wind tunnel testing will 
validate this result. 

3. Grid Studies 

A grid density study was performed on the baseline axisymmetric Nozzle 6 configuration. The grid study 
reduced the number of gridpoints for all zones by a factor of 2 and then a factor of 3. The reduction of gridpoints by 
a factor of 2 had a measurable effect on the AP/Poo signature at 150 in. from the nozzle centerline, reducing the 
magnitude of peak pressures by 10 percent (see Fig. 15). Further reduction to a factor of 3 was similar. Grid density 
was then increased by a factor of 1 .4, and results were nearly identical to the original baseline grid. As a result, all 
CFD cases kept the original dense grid at 51 1,299 gridpoints. 

4. Nozzle Boat Tail Angle 

The nozzle boat tail angle had an effect on the initial expansion shock wave from the nozzle. To study this effect, 
the baseline nozzle was modified to create a 2-degree boat tail angle. Figure 16 makes a direct comparison of nozzle 
boat tail angle on the boom signature for an NPR of 8, and shows the reduction in boom signature. 

For the baseline Nozzle 6, results for near-field pressures at 150 in. from the nozzle centerline are summarized in 
Fig. 17. Results are for NPRs from 7 to 9, and are more difficult to interpret, as the signatures are hard to 
distinguish. To clarify differences, the estimated boom signature was computed and is shown in Fig. 18. Differences 
are noticeable between an NPR of 7 and 8. For NPR of 9, the boom signature becomes stronger. For the 2-degree 
boat tail angle, there is a narrower band of pressure ratios, between 7 and 8, where there is a benefit in adjusting the 
nozzle pressure ratio. Further study can determine an optimum configuration of boat tail angle and nozzle operating 
parameters. 




400 450 500 550 

x 

Fig. 15. Baseline Nozzle 6 grid study results. 
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Fig. 16. Estimated sonic boom signature at ground through an ideal 
atmosphere, computed with PCBoom. Nozzle 6 comparison of 5-degree and 
modified 2-degree boat tail angle at NPR 8. 



Fig. 17. Far-field pressure profiles AP/Poo at y = 150 in. from nozzle 
centerline. Nozzle 6 with a modified 2-degree boat tail angle at NPR 7 to 9. 
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Fig. 18. Estimated sonic boom signature at ground through an ideal 
atmosphere, computed with PCBoom. Nozzle 6 with a modified 2-degree 
boat tail angle at NPR 7 to 9. 



NPR 


Fig. 19. Baseline Nozzle 6 performance curve. 


III. Conclusions 

CFD analysis was used to address the plume effects from isolated exhaust nozzles on sonic boom signature. The 
CFD analysis used Nozzle 6 from NASA TN-D-5553 as a baseline. Values of AP/P^ were computed for x/D values 
less than 1 0 and at 1 nozzle diameter from the nozzle centerline, and agreement between the CFD analysis and the 
NASA TN-D-5553 experimental data for Nozzle 6 was acceptable. Results show that there was a second AP/P 
overpressure beyond an x/D of 10 that was not captured in NASA TN-D-5553. 

The API Poo pressure signatures at 150 in. from the nozzle centerline were collected for the baseline axisymmetric 
Nozzle 6 configuration. Pressure signatures were collected for nozzle pressure ratios between 6 and 10. PCBoom 
extrapolated these signatures to a ground-observed sonic boom N-wave. Trends show a reduction in sonic boom N- 
wave signature as NPR is increased from 6 to 10, as underexpanded nozzle operation reduces the boat tail expansion 
and recompression shock. Nozzle boat tail angle has an important effect on the reduction of sonic boom signature; 
further study is needed to optimize combinations of nozzle pressure ratio and boat tail angle. For the pressure ratios 
studied, the performance curve for this nozzle is flat (see Fig. 19), so there is no significant loss in thrust coefficient 
as the NPR is increased. As a result, this benefit could be realized without significant loss of performance. 
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The AP/Pw pressure signatures at 150 in. from the nozzle centerline were also collected for a high aspect ratio 
nozzle based on the Nozzle 6 design for comparison. Pressure signatures were collected for nozzle pressure ratios 
from 8 to 12. The AP/P X signatures were nearly twice as strong for the two-dimensional case, and trends also show a 
reduction in sonic boom signature as NPR is increased from 8 to 12. 

As future low boom aircraft designs are developed and improved, there will be a need for improved 
understanding of the interaction between the aircraft boat tail shocks and the exhaust nozzle plume. These CFD 
analyses are intended to provide a baseline study for future analysis efforts. Future efforts will be needed to 
investigate plume interaction trends on an integrated aircraft and propulsion system problem. 
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